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We study the effect of quantum memory in non-inertial frames under the influence of 
amplitude damping, depolarizing, phase flip and bit-phase flip channels. It is shown that 
«) ' the entanglement of initial state is heavily influenced by quantum correlations. It is seen 

that quantum memory compensates the loss of entanglement caused by the Unruh effect. It 
i TO . is interesting to note that the sudden death of entanglement disappears for any acceleration 



in 



for higher values of quantum memory. Therefore, it is possible to avoid ESD in non-inertial 
frames due to the presence of quantum memory. Furthermore, the degree of entanglement 



ph \ is enhanced as we increase the degree of memory and it maximizes for maximum correlations. 

^ ' Keywords: Quantum memory; entanglement; quantum channels. 
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I. INTRODUCTION 



Quantum entanglement being a fundamental potential resource for communication is one of the 
key quantitative notions of the intriguing field of quantum information and quantum computation. 
Entanglement behavior in non-inertial frames was first considered by Alsing et al. [1]. They studied 
the fidelity of teleportation between relative accelerated partners. Afterwards, a number of authors 
have focused on the topic of the quantum information in a relativistic setup [2-8]. Schuller et al. 
[9] have studied how the Unruh effect changes the degree of quantum entanglement. However, 
all the investigations in non-inertial frames were confined to the studies of quantum information 
in an isolated system except [10, 11], where the authors have analyzed the decoherence effects in 
non-inertial frames. 

Since in a realistic quantum system, the interaction between the quantum system and the 
surrounding environment is inevitable, it is therefore, important to study the effect of environmental 
influences on the entanglement dynamics of a system. Quantum channels can be implemented by 
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suitable quantum devices consisting of intrinsic degrees of freedom associated with the environment 
and acting on the system via particular interactions between the system and the environment. The 
assumption that noise is uncorrelated between successive uses of a channel is not realistic. Hence 
memory effects need to be taken into account. Quantum channels with memory [12-14] are the 
natural theoretical framework for the study of any noisy quantum communication system where 
correlation time is longer than the time between consecutive uses of the channel. A more general 
model of a quantum channel with memory was introduced by Bowen and Mancini [15] and also 
studied by Kretschmann and Werner [16]. 

In this paper, we investigate the effect of quantum memory on the entanglement of Dirac fields 
in non-inertial frames. We consider the amplitude damping, depolarizing, phase flip and bit-phase 
flip channels. We assume that two observers, Alice and Rob, share an entangled initial state at 
the same point in flat Minkowski spacetime. After that Alice stays stationary while Rob moves 
with uniform acceleration. We assume that Alice and Rob have detectors that are sensitive only 
to their respective modes and both the observers share the following maximally entangled initial 
state 

|*>Afl = ^(|0> J 4|0)fl+|l)A|l>Jl) (1) 

where the two modes of Minkowski spacetime that correspond to Alice and Rob are \t})a and \rj)n 
respectively. Prom the perspective of Rob, the Minkowski vacuum is found to be a two-mode 
squeezed state [4] 

|0) M = cosr|0)j|0)jj + sinr|l)i|l)jj (2) 

where cosr = ( e - 27I W a + ly 1 / 2 . The constants w, c and a, in the exponential stand for Dirac 
particle's frequency, light's speed in vacuum and Rob's acceleration respectively. The subscripts 
/ and II of the kets represent the Rindler modes in region I and II, respectively, in the Rindler 
spacetime diagram (see Ref. [10], Fig. (1)). The excited state in Minkowski spacetime is related 
to Rindler modes as follows [4] 

|1)m = |1)/|0)// (3) 

Using equations (2) and (3), equation (1) can be written in terms of Minkowski modes for Alice 
and Rindler modes for Rob as 

\V) A ,I,II = -4 (coBr|0) A |0)/|0)//| + sinr|0) A |l)/|l)// + |1)a|1)j|0>//) (4) 
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Since Rob is causally disconnected from region II, therefore, by tracing over all the modes in region 
II leaves the following mixed density matrix between Alice and Rob 



We study the entanglement of Dirac fields in non-inertial frames influenced by different memory 
channels, such as amplitude damping, depolarizing, phase flip and bit-phase-flip channels, param- 
eterized by decoherence parameter p and memory parameter p. Here p £ [0, 1] and p £ [0, 1] 
represent the lower and upper limits of decoherence parameter and memory parameter respec- 
tively. It is seen that depolarizing channel influences the entanglement of Dirac fields heavily in 
comparison to the amplitude damping and flipping channels. Therefore, the entanglement of Dirac 
fields is strongly dependent on degree of correlations of the noisy channels. 



Since noise is a major hurdle while transmitting quantum information from one party to other 
through classical and quantum channels. This noise causes a distortion of the information sent 
through the channel. If multiple uses of a channel are not correlated, there is no advantage in using 
entangled states. Correlated noise, also referred as memory in the literature, acts on consecutive 
uses of the channel. However, in general, one may want to encode classical data into entangled 
strings or consecutive uses of the channel may be correlated to each other. Hence, we are dealing 
with a strongly correlated quantum system, the correlation of which results from the memory of 
the channel itself. In Ref. [12], a Pauli channel with partial memory was studied. The action of 
such a channel on two consecutive qubits is given in Kraus operator form as 



where Gi (crj) are usual Pauli matrices, pi (pj) represent the quantum noise and indices i and j runs 
from to 3. The above expression means that with probability \x the channel acts on the second 
qubit with the same error operator as on the first qubit, and with probability (1 — \i) it acts on 
the second qubit independently. Physically the parameter p is determined by the relaxation time 
of the channel when a qubit passes through it. In order to remove correlations, one can wait until 
the channel has relaxed to its original state before sending the next qubit, however this lowers the 
rate of information transfer. Thus it is necessary to consider the performance of the channel for 
arbitrary values of /i to reach a compromise between various factors which determine the final rate 



p AI = I (cos 2 r|00) (00| +cosr(|00) (11| + |11) (00|) + sin 2 r|01) (01| + |11) (11|) (5) 



II. ENTANGLEMENT IN A CORRELATED ENVIRONMENT 




(6) 
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of information transfer. Thus in passing through the channel any two consecutive qubits undergo 
random independent (uncorrelated) errors with probability (1 — fx) and identical (correlated) errors 
with probability /i. This should be the case if the channel has a memory depending on its relaxation 
time and if we stream the qubits through it. A detailed list of single qubit Kraus operators for 
different quantum channels with uncorrelated noise is given in table 1. The action of such a channel 

if n qubits are streamed through it, can be described in operator sum representation as [17] 

n-i 

pf= E a® ^xkn^L® 4j (7) 

fcl i ....,.fcn = 

where p in represents the initial density matrix for quantum state and Ak n are the Kraus operators 
expressed in equation (6). The Kraus operators satisfy the completeness relation 

E 4a. = 1 (8) 



However, the Kraus operators for a quantum amplitude damping channel with correlated noise are 
given by Yeo and Skeen [13] as given as 
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00 



A c 



cos x 
10 
10 

sinx 

where, < x ^ 7r /2 and is related to the quantum noise parameter as 

sin x = VP 









(9) 



(10) 



It is clear that Aq cannot be written as a tensor product of two two-by-two matrices. This gives 
rise to the typical spooky action of the channel: 1 01) and |10), and any linear combination of them, 
and 1 11) will go through the channel undisturbed, but not 1 00) . The action of this non-unital 
channel is given by 



7T 



P 



3>(tt) = (1 - M ) A l^ U ii + V- E A kk«A 

i,j=0 



(11) 



k=0 



The action of the super-operators provides a way of describing the evolution of quantum states 
in a noisy environment. In our scheme, the Kraus operators are of the dimension 2 2 . They 
are constructed from single qubit Kraus operators by taking their tensor product over all n 2 
combinations 



Ai 



(12) 
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where n is the number of Kraus operator for a single qubit channel. The final state of the system 
after the action of the channel can be obtained as 



Pf = ®pApa,i) 



(13) 



where is the super-operator realizing the quantum channel parametrized by real numbers p 
and fj,. The density matrix after the action of different memory channels are obtained by using 
equations (5-13). It is important to mention here that we consider the case that both the qubits 
are coupled to the time correlated channel with memory. The final density matrix after the action 
of amplitude damping channel is given by 

/ !(-(-! + n)p(l+p) 



n AD 
Pf 



-(— 1 + p) cos 2 r) 




!(i_ p+Ai (_i 





i(l + (-l + ^ 2 
+(—1 + p — \ip) cos 2 r) 









5(-l + M) 
(-l+p)p 





1(1_P +M (_1 



\ 



^ - p) +p))cosr 



+ a/(1 — p) +p))cosr 




i (i_(_i + M ) 

(— 2 + p)p + /ipcos 2 r) , 
(14) 



The final density matrix after the action of depolarizing channel is given by 



Dcp 

Pf 



|cos 2 r(2 + p(-3 
-2/i(-2 + p) + 2p) 
-(-2+p)(l + (-l 
+/x)p) cos 2 r) 





| cos r(2 + p(— 3 
-2/x(-2 + p) + 2p) 
-(-2+p)(l + (-l 
+/i)p) cos 2 r) 



^(4 + (-2 + M )p 
+2(-2 + p)(l + (-l 
+/i)p) cos 2 r) sin 2 r 




-jQ^pcosr sin 2 r 



^cosr(4+(-2 \ 

+/x)p + (4+p(-10 

+/i(H - 6p) +6p)) 
cos 2 r) 



±(4 + (-2 + f i)p+ 
(4 + p(-10 + /i(ll 
— 6p) + 6p)) cos 2 r) 



(15) 



The final density matrix after the action of bit-phase flip channel is given by 



Bpf 
Pf 



( i(l+ p (-l-2/i(-l+p) 
+2p)-(-l+p)(l 
+2(-l +p)p) cos(2r)) 









| cos r(2 + p(— 3 — 2p 
(-2 + p) + 2p) - (-2+p) 
(1 + (-1 + m)p) cos 2 r) 



^cosr(4+ (-2 + p)p 
+(4 + p(-10 + /i(ll 
— 6p) + 6p)) cos 2 r) 



^(4 + (-2 + /i)p+ 
2(-2+p)(l + (-l 
cos 2 r) sin 2 r 




-is/*? 

cos r sin 2 r 





^(4+(-2 + /i)p 
+(4 + p(-10 + /i(ll 
— 6p) + 6p)) cos 2 r) 



(16) 



The final density matrix after the action of phase flip channel is given by 



Pf 
Pf 



1 2 

2 cos r 



^ sin 2 r —jQ^ip cos r sin 2 r 



|(l + 4p(-l+ 
V p + p — //p)) cos r 













i(l + 4p(-l + /x \ 
+p — pp)) cos r 







(17) 



/ 



where the super-scripts AD, Dep, Bpf and Pf correspond to the amplitude damping, depolarizing, 
bit-phase flip and phase flip channels respectively. In order to study the degree of entanglement 
in the two qubits mixed state in a correlated noisy environment we calculate the concurrence C, 
which is given as [18] 

(18) 



C = max{0, y/Tx - ^f\~ 2 - \/A3 - y%} A* ^ A;+i ^ 



where Aj are the eigenvalues of the matrix pffif, with pj being the spin flip matrix of pj (as given 
in equations (14-17)) and is given as 



Pf = i a y ® a y)Pf( a y ® a y) 



(19) 



7 



where a y is the usual Pauli matrix. The eigen values of pj pj are given as under 
1 



2 2 2 2 3 3 2 3 1 

-(p — [ip — p + 3 fip — 2{ip — p + 2 fip -/i p + p 
4 

— 2\xp 4 + (i 2 p 4 + 2 cos 2 r — 2fi cos 2 r + 2fi 2 cos 2 r + 2^-^/1 — p cos 2 r 
— 2fj? — p cos 2 r — bp cos 2 r + 6/ip cos 2 r — 3(i 2 p cos 2 r — 2\i\J 1 — pp cos 2 r 
+2// 2 yl—pp cos 2 r + 4p 2 cos 2 r — 4:fip 2 cos 2 r — p 3 cos 2 r + 2//p 3 cos 2 r 
—fi 2 p 3 cos 2 r + fip cos r 4 — pp 2 cos r 4 



V2 



^3,4 



((-1 +p)(-l +p- 2/x(-l + v / T^P + p) 
+/x 2 (-2 + 2 VI -P + p)) cos 2 r 
(4 - 4p + + 4p 2 + 13 fip 2 - 20/ii V 
-12p 3 + 24//p 3 - 12/i 2 p 3 + 8p 4 - 16 fip 4 + 8fi 2 p 4 
-4(-l - 3(-l + fi)p + (-3 + 3/i + ^ 2 )p 2 
\ +(-l + At)V)cos(2r) - +p)pcos(4r))))) 



-(-1 + +P)p(l + (-1 + /")P 2 + (-1+p-MP) cos2r ) 



Similarly, the eigen values of p^ ep p^ ep are given below 

1 



Ai = — cos 2 r(8 + 2(-8 + 9/i)p +(14-19// + 4fi 2 )p 2 
32 

-2(2 - 3fi + ^ 2 )p 3 + (16 + 48(-l + n)p 
+2(30 - 52/z + 23/i 2 )p 2 + (-40 + 87/i - 47/i 2 )j> 3 
+12(-1 + /x) 2 p 4 ) cos 2 r + (-2 + p)(-4 + (14 - 15/i)p 
+ (-16 + 27 fi - H/i 2 )p 2 + 6(-l + //) V) cosr 4 ) 




A 



(20) 



2,3,4 



(21) 
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and the eigen values of Pj Pf Pj Pf are obtained as 



Ai,: 



1 

32 ( 



27p + 30/ip + 55p 2 - 86/ip 2 + 28/i V - 56p 3 



+112//p 3 - 56/iV + 28p 4 - 56/V + 28/i 2 p 4 

+8cos(2r) - 36pcos(2r) + 32/zpcos(2r) + 68p 2 cos(2r) 

-96//p 2 cos(2r) + 32/x V cos(2r) - 64p 3 cos(2r) 

+128/V cos(2r) - 64/i 2 p 3 cos(2r) + 32p 4 cos(2r) 

— 64^tp 4 cos(2r) + 32/i 2 p 4 cos(2r) — pcos(4r) 

+2/ipcos(4r) + 5p 2 cos(4r) — 10/ip 2 cos(4r) 

+4/i 2 p 2 cos(4r) - 8p 3 cos(4r) + 16/ip 3 cos(4r) 

— 8/i 2 p 3 cos(4r) + 4p 4 cos(4r) — 8^p 4 cos(4r) + 4// 2 p 4 cos(4r) 



±4\/2 



((1 + 2(-l + /x)p - 2(-l + ^p 2 ) 2 cos 2 r(4 - lip 
+14//p + 23p 2 - 38/ip 2 + 12/i V 
-24p 3 + 48/xp 3 - 24//V + 12p 4 - 24^p 4 
+12//V + 4(1 + (-5 + 4/x)p + (3 - 2fi) V 
-8(-l + ^p 3 + 4(-l + ju) V) cos(2r) 
\j +(-1 +p)p((l - 2p) 2 - 2//(l - 2p) 2 + 4/i 2 (-l + p)p) cos(4r)))) 



1 



2^3 



A 3j4 = — (5p - 2/xp + 23p 2 - 54/ip 2 + 28/*V - 56p 3 + 112/zp d - 56/Tp 

+28p 4 - 56pp 4 + 28// 2 p 4 - 4pcos(2r) + 36p 2 cos(2r) - 64^p 2 cos(2r) 
+32/i 2 p 2 cos(2r) - 64p 3 cos(2r) + 128pp 3 cos(2r) - 64/x 2 p 3 cos(2r) 
+32p 4 cos(2r) — 64/ip 4 cos(2r) + 32p 2 p 4 cos(2r) — pcos(4r) 
+2//pcos(4r) + 5p 2 cos(4r) — 10/ip 2 cos(4r) + Afi 2 p 2 cos(4r) 
-8p 3 cos(4r) + 16/ip 3 cos(4r) - 8fi 2 p 3 cos(4r) 
+4p 4 cos(4r) — 8[ip 4 cos(4r) + 4p 2 p cos(4r) 



±8^ 



((-1 + /i) 2 (-l + p) 3 p 3 cos 2 r(-5 + 2/i - 12p + 24/xp 
-12/i 2 p + 12p 2 - 24pp 2 + 12/i 2 p 2 
+4(1 - 4(-l + p) 2 p + 4(-l + (i) V) cos(2r) 
\J +((1 - 2p) 2 - 2/i(l - 2p) 2 + 4 / u 2 (-l + p)p) cos(4r)))) 



(22) 
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The eigen values of 



PfPf are given as under 



A 2 



Ai 



(1 + 2(-l + fi)p - 2(-l + p)p 2 ) 2 cos 2 r 
4(-l + p) 2 (-l +p) 2 p 2 cos 2 r 








(23) 



The concurrence C, is calculated using equation (18). It can be easily checked from equations 
(20-23) that the concurrence becomes cosr if we set p = p = 0, which reproduces the results of 
Ref. [4]. Similarly the results of Ref [10-11] can be reproduced by setting p = 0. 



We have computed the mathematical relations for concurrence in case of amplitude damping, 
depolarizing, phase flip and bit-phase flip channels. However, due to very lengthy mathematical 
expressions, we have only plotted it for different parameters such as p, r, p respectively. 

In figures 1 and 2, we plot the concurrence as a function of memory parameter p for different 
values of Rob's acceleration, r for decoherence parameter p = 0.5 for amplitude damping, depolar- 
izing channels respectively. In figures 3 and 4, we plot the concurrence as a function of memory 
parameter p with others parameters similar as in figures 1 and 2 for bit-phase flip and phase flip 
channels respectively. It is seen that quantum memory enhances the entanglement of Dirac fields 
in contrast to the effect of decoherence. It is also seen that the depolarizing channel influences 
the entanglement more heavily in comparison to the other channels. The depolarizing channel 
enhances the degree of entanglement upto 1 even with 50% decoherence (see figure 2). It is clear 
from figure 3 that for lower values of memory parameter, the decoherence dominates and ESD 
behaviour can be seen. However for higher values of memory, the ESD can also be avoided for 
flipping channels (see figures 3 and 6). 

In figure 5, we plot the concurrence as a function of Rob's acceleration r for p = p = 0.5 
for amplitude damping, depolarizing, bit-phase flip and phase flip channels. It is seen that the 
entanglement degradation is controlled due to the presence of quantum memory even at 50% 
decoherence. It is also seen that entanglement is heavily degraded by amplitude damping channel 
as compared to other channels. In figure 6, we plot the concurrence as a function of decoherence 
parameter p for p = 0.5 and r = ir/6 for amplitude damping, depolarizing, bit-phase flip and 
phase flip channels. It is seen that the maximum entanglement degradation occurs for amplitude 
damping channel as compared to other channels. Furthermore, the symmetric behaviour of flipping 



III. DISCUSSIONS 



10 



channels is seen with maximum degradation for 50% decoherence. In order to check the validity of 
our calculations, we plot in figure 7, the concurrence as a function of decoherence parameter p for 
[i = and r = 7r/10 for amplitude damping, depolarizing, bit-phase flip and phase flip channels. It 
is seen that the maximum entanglement degradation occurs for bit-phase flip channel as compared 
to other channels. It is clear from the figure that the graphs of Ref. [10, 11] are reproduced in the 
absence of memory. 



IV. CONCLUSIONS 



We study the effect of quantum memory in non-inertial frames for amplitude damping, 
depolarizing, phase flip and bit-phase flip channels. It is shown that the entanglement of initial 
state is heavily influenced by quantum correlations. It is seen that quantum memory compensates 
the loss of entanglement generated by the Unruh effect and decoherence. It is interesting to 
note that the sudden death of entanglement disappears for any acceleration for higher degree of 
memory. Furthermore, the entanglement of Dirac fields increases as we increase the degree of 
memory. In conclusion, we can say that quantum memory can avoid ESD caused by Unruh effect 
for all the channels considered. 
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Figures captions 

Figure 1. The concurrence is plotted as a function of memory parameter fi for different values 
of Rob's acceleration, r = (solid line), r = ir/6 (dashed line) and r = tt/4 (dotted line) for 
decoherence parameter p = 0.5 for amplitude damping. 

Figure 2. The concurrence is plotted as a function of memory parameter fi for different values 
of Rob's acceleration, r = (solid line), r = ir/6 (dashed line) and r = tt/4 (dotted line) for 
decoherence parameter p = 0.5 for depolarizing channel. 

Figure 3. The concurrence is plotted as a function of memory parameter fi for different values 
of Rob's acceleration, r = (solid line), r = 7r/6 (dashed line) and r = 7r/4 (dotted line) for 
decoherence parameter p = 0.5 for bit-phase flip channel. 

Figure 4. The concurrence is plotted as a function of memory parameter \x for different values 
of Rob's acceleration, r = (solid line), r = ir/6 (dashed line) and r = ir/4 (dotted line) for 
decoherence parameter p = 0.5 for phase flip channel. 

Figure 5. The concurrence is plotted as a function of Rob's acceleration r for p = [i = 0.5 for 
amplitude damping (solid line), depolarizing (dashed line), bit-phase flip (dotted line) and phase 
flip (dot dashed line) channels. 

Figure 6. The concurrence is plotted as a function of decoherence parameter p for fi = 0.5 and 
r = 7r/6 for amplitude damping (solid line), depolarizing (dashed line), bit-phase flip (dotted line) 
and phase flip (dot dashed line) channels. 

Figure 7. The concurrence is plotted as a function of decoherence parameter p for /x = and 
r = 7r/10 for amplitude damping (solid line), depolarizing (dashed line), bit-phase flip (dotted 
line) and phase flip (dot dashed line) channels. 

Table Caption 

Table 1. Single qubit Kraus operators for amplitude damping, depolarizing, bit-phase flip and 
phase flip channels where p represents the decoherence parameter. 
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FIG. 1: The concurrence is plotted as a function of memory parameter /j, for different values of Rob's 
acceleration, r = (solid line), r = ir/6 (dashed line) and r = tt/4 (dotted line) for decoherence parameter 
p = 0.5 for amplitude damping. 
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FIG. 2: The concurrence is plotted as a function of memory parameter /j, for different values of Rob's 
acceleration, r = (solid line), r = ir/6 (dashed line) and r = tt/4 (dotted line) for decoherence parameter 
p = 0.5 for depolarizing channel. 
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FIG. 3: The concurrence is plotted as a function of memory parameter \x for different values of Rob's 
acceleration, r = (solid line), r = ir/6 (dashed line) and r = tt/4 (dotted line) for decoherence parameter 
p = 0.5 for bit-phase flip channel. 
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FIG. 4: The concurrence is plotted as a function of memory parameter /j, for different values of Rob's 
acceleration, r = (solid line), r = ir/6 (dashed line) and r = tt/4 (dotted line) for decoherence parameter 
p = 0.5 for phase flip channel. 
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FIG. 5: The concurrence is plotted as a function of Rob's acceleration r for p = = 0.5 for amplitude 
damping (solid line), depolarizing (dashed line), bit-phase flip (dotted line) and phase flip (dot dashed line) 
channels. 
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FIG. 6: The concurrence is plotted as a function of decoherence parameter p for (i = 0.5 and r = tt/6 for 
amplitude damping (solid line), depolarizing (dashed line), bit-phase flip (dotted line) and phase flip (dot 
dashed line) channels. 
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FIG. 7: The concurrence is plotted as a function of decoherence parameter p for /i = and r = tt/10 for 
amplitude damping (solid line), depolarizing (dashed line), bit-phase flip (dotted line) and phase flip (dot 
dashed line) channels. 
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TABLE I: Single qubit Kraus operators for amplitude damping, depolarizing, bit-phase flip and phase flip 
channels where p represents the decoherence parameter. 



Amplitude damping channel 


A = 


I 

y/l-p 


,A 1 = 


o Vp 




Depolarizing channel 


A = ^/l-%I, A 1 = ^a x 


Bit-phase flip channel 


Aq = y/1 - pi, Al = y/p(Ty 


Phase flip channel 


A = V 1 - Pi, A i = Vp^z 



